Quasi-neutral limit of the Euler-Poisson and 
Euler-Monge- Ampere systems 



G. Loeper* 



Abstract 

This paper studies the pressureless Euler-Poisson system and its fully non-linear counter- 
part, the Euler-Monge- Ampere system, where the fully non- linear Monge- Ampere equation sub- 
stitutes for the linear Poisson equation. While the first is a model of plasma physics, the second 
is derived as a geometric approximation to the Euler incompressible equations. Using energy 
estimates, convergence of both systems to the Euler incompressible equations is proved. 

1 Introduction 

In this paper we consider a model of a collisionless plasma where the ions are supposed to be at 
rest and create a neutralizing background field. The motion of the electrons can then be described 
by using either the kinetic formalism or the hydrodynamic equations of conservation of mass and 
momentum as we do here. The self-induced electric field is the gradient of a potential that depends 

on the electron's density p either through the linear Poisson equation: A0 = — (p — 1), or through the 

e 

fully non-linear Monge- Ampere equation: det(J + edij<p) = p. This gives the Euler-Poisson ((EP e )) 
system and Euler-Monge- Ampere ((EMA e )) system. The non-dimensional rescaled version of both 
systems is the following: 

dtp + V • (pv) = 0, 

d t v + v ■ Vt> = -V0, 
e 

eA0 = p — 1 in the Poisson case , 

det(J + edijtfi) = p in the Monge- Ampere case. 

Note that the systems are pressureless, and the only force is due to electrostatic interaction. The 
energy of those systems is given by 



£€ = IJ Pl^| 2 + |V0| 2 dx 



"Research supported by a doctoral grant at Laboratoire J. A. Dieudonne, Universite de Nice-Sophia- Antipolis 



G. Loeper, Euler-Poisson and Euler-Monge- Ampere 



2 



for Euler-Poisson and 

£t = \j pM 2 + pIWI 2 dx 

for Euler-Monge-Ampere. 

The asymptotic we look at consists in considering large scales compared to the Debye length (e). 
At those scales the plasma appears to be electrically neutral. In this limit the plasma is expected 
to behave like an incompressible fluid, therefore governed by the incompressible Euler equation (E). 
We intend to rigorously justify those limits in the present work. 

Physical interpretation of the quasi-neutral limit for Euler-Poisson The complete model 
of collisionless plasma describes the behavior of two species: the ions and the electrons. However the 
ratio of the electron's mass and the ion's mass is of several orders of magnitude, therefore we make 
the assumption that the ions are at rest, and distributed over a regular grid. This assumption will 
imply the neutralizing background electric field (the -1' term in the Poisson equation eA0 = p — 1). 
The parameter e comes from the vacuum permittivity, obtained after many rescalings of the equation. 
The typical value of e 2 is between 10~ 10 and 10~ 5 . We consider therefore e as a small parameter, 
and investigate the limit e — > of the Euler-Poisson system. Note that if (v, p, <p) is a solution to the 
Euler-Poisson system with e = 1, (v t ,p t ,4> t ) := (v , p, (p)(e~H, e _1 x) is a solution of the Euler-Poisson 
system with parameter e. Therefore, the limit e — ► can be interpreted as a study of the long time 
- large scale behavior of the system. 

Geometric interpretation of the quasi-neutral limit for Euler-Monge-Ampere Whereas 
the Euler-Poisson system relies on a well known physical model, the Euler-Monge-Ampere system, 
less famous, is a fully non-linear (but asymptotically close in the quasi-neutral regime) version of 
the Euler-Poisson system; it can be seen as a non-linear model of electrostatic interaction with the 
advantage of allowing finite electric field for point charges (see also [B] where the Born-Infeld system 
of electromagnetism is studied, a system that exhibits similar non-linear features). Apart from 
this interpretation, the main motivation for the study of the Euler-Monge-Ampere system is the 
following: it appears as a 'canonical' relaxation of the geodesies on the group of measure preserving 
diffeomorphisms (therefore of the Euler incompressible equation, see [2]). This interpretation will be 
developed more accurately in the section 13.11 devoted to the Euler-Monge-Ampere equation. This 
model was first introduced, in a discrete version, by Y. Brenier in j^j. Later, a kinetic version, the 
Vlasov-Monge- Ampere system, has been introduced and studied by Y. Brenier and the author in [Zj. 
The present work may be seen as a further step in this study. 

To see why the (EP e ) and the (EMA e ) systems should be asymptotically close in the quasi-neutral 
regime, notice that if p is close to 1 then should be small, hence det (I + edij (ft) = 1 + eA0 + O(e 2 ) 
and one recovers the Poisson equation. For this reason the proof of the convergence of both systems 
will be very close and this is why we present them altogether. 
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Related results concerning singular perturbations This work is concerned with the motion 
of slightly compressible fluids seen as singular perturbations of the Euler incompressible equation; 
this field has been widely investigated using different techniques: 

- Traditional analysis and geometry on the group of volume preserving diffeomorphisms (see 
section l3~T|) as done in ^Hj where the convergence holds in H s norm, s large, for well prepared initial 
data, restricted to the case of barotropic fluids (i.e. when the pressure is a local function of the 
density p, a case different from the one studied here.) 

- Energy estimates as done in pH] again for the case of barotropic fluids, where convergence holds 
in all H s norms for well prepared initial data. The result has also been extended to the non-isentropic 
case by |2T| . 

- Pseudo-differential energy estimates as done in [16J which can be seen as a pseudo-differential 
generalization of JB] and where the same convergence holds for a broader class of singular perturba- 
tions, including non-local dependence between density and pressure. 

- Modulated energy techniques for convergence of the Vlasov-Poisson system to the so-called 
dissipative solutions of the Euler equation, as done in ^] , and for the Vlasov-Monge- Ampere system 
(a kinetic version of the Euler-Monge- Ampere system) as done in jjj. The convergence result obtained 
there holds in weighted "L 2 " norms, this method has the advantage to be valid for weak solutions, 
and does not require any smoothness of the solution. We give here more details on this result: 

Modulated energy technique for Vlasov-Poisson and Vlasov-Monge-Ampere Those sys- 
tems are the kinetic extensions of the (EP e ) and (EMA e ) systems. They read as follows: 



In ^] and [7j, the following results have been obtained: 

Theorem 1.1 Let f be a weak solution 0/ J3 0) (resp. 0/ H ||)J with finite energy, let 

(t,x) — > v(t,x) be a smooth solution of the incompressible Euler equation (0) for t G [0,T], and 
p(t, x) the corresponding pressure, let 



where E e (t) = e 2 J p|V0| 2 /2 (resp. E e (t) = e 2 J |V0| 2 /2 in the Poisson case), then 

H e (t) < C exp(Ct)(H e (0) + Ce 2 ), Vt G [0,T]. 
The constant C depends only of the W^'°° norm of {v(s, -),p(s, .),d t p(s, .),'Vp(s, .) s G [0,T]}. 



(1) 




(2) 

(3) 
(4) 




e 
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Results Here we shall obtain by energy estimates a convergence to the Euler incompressible sys- 
tem in L^H S X norm for any s large enough (the minimal smoothness will be made precise). The 
convergence of both systems holds on the range of time on which the solution of Euler is smooth 
enough (roughly speaking, we will need at least D 2 v to be bounded in L°°). Our work is based on 
the modulated energy techniques, restricted to the case of monokinetic velocity profiles. Indeed, the 
quasi-neutral limit is much more difficult without this assumption, and is even known to be false in 
some cases (the two-streams instability, see [T2j). We will mostly restrict ourselves to the case of 
well-prepared initial data, but we will investigate briefly the case of non-prepared initial data: in that 
case the divergence part of the initial velocity is not assumed to be small, and we only assume that 
the initial fluctuations of the electronic density p are of order e, so that the energy remains bounded. 
The electric field is expected to oscillate at frequency e _1 and with amplitude 0(e _1 ). This case will 
be treated in section 2.3, and we will obtain that the divergence-free part of the velocity converges 
strongly to a solution of the Euler incompressible equation while its potential part stays bounded, 
but oscillates strongly with respect to time, and therefore converges weakly to 0. 

We obtain also that both systems are closer to each other than they are close to the Euler 
incompressible system; (EP e ) is thus a corrector in the convergence of (EMA e ) to (E). 

Although the operators that define the acceleration from the density are differential operators 
(and even fully non-linear in the second case), our proof does not use the pseudo-differential formal- 
ism. Actually, we were not able to use the general theorem obtained by Grenier pB] for singular 
perturbations: it might be because of the absence of pressure that changes the symmetrizers of the 
system. After a convenient change of variable however, the system appears under a form which is 
strongly reminiscent (at least for the highest order terms) of the rapidly rotating fluids. This limit 
has been treated in [To] . 

Finally we also mention the work of Cordier & Grenier [IT], and Wang where the quasi- 
neutral limit 'with pressure' is treated. The techniques used there do not apply here, and it is worth 
noting that the scaling obtained are not the same. The reader can also refer to ^2] and j22] where 
different regimes of the Euler-Poisson system are studied. 

We split the rest of the paper in two sections: the first one devoted to the study of the Euler- 
Poisson system, and the second devoted to the study of the Euler-Monge-Ampere system. 

1.1 Notation 

Hereafter x G T d = R d /Z d and t G R + ; v(t,x) G R d stands for the velocity and p(t,x) G R + is the 
macroscopic density of electrons; <fi(t,x) G R is the electrostatic potential; d = 2 or 3. 
It is always assumed that pit, ■) has total mass equal to 1. 

The divergence of a vector field v will be denoted by divt> or V • v; its rotational (or curl) will be 
denoted by V x v or curlf . 

The components of a vector will be denoted with superindices, i.e. v G R d = (v l )i=i..d- 
In all the paper, [•] will denote the integer part. 

We denote respectively by (E), (EP e ), (EMA e ) the systems Euler incompressible, Euler-Poisson, 
Euler-Monge-Ampere. 
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2 The Euler-Poisson system 

We consider the following Euler-Poisson system denoted by (EP e ): 

(5) d t p + V ■ (H = 0, 

(6) d t v + v • Vv = ^— , 

(7) eA<f> = p-l, 

and consider the limit e going to 0. We recall also the incompressible Euler equation (E): 

d t v + v ■ Vf = Vp, 

(8) V • v = 0. 

We recall (see jH] for a reference on the topic) that in the periodic case, the Cauchy problem for (JHJ) is 
well posed in H s (T d ) if s > d/2 + 1. More precisely, if d = 2, for any divergence-free initial datum in 
H s , there exists a unique global solution in L^ C (R, if s (T d )); if d = 3 one can only prove the existence 
of a smooth solution in finite time, belonging to L^ c ([0, T[, H s (T d )) for some T > 0. We will then 
prove the following: 

2.1 Result 

Theorem 2.1 Let sGN wift s > [d/2]+2. Letvo be a divergence- free vector field on T d . Let (v,p) be 
a smooth solution of the Euler incompressible system |2J) on [0, T] x T d , with initial data vq, satisfying 
v E L°°([0,T], H s+1 (T d )) . Let (Pq,Vq) be a sequence of initial data such that J Td p e (x)dx = 1, and 
such that (e -1 (i>o — v ), e~ 2 (p e — 1)) is bounded in H s x H s ~ 1 (Y d ). Then there exists a sequence 
(v e ,p e ) of solutions to (EP e ) with initial data (vq,p € ), belonging to L°°([0, T e ], H s x H^ 1 ^)), with 
liminf^o^'e > T. Moreover for anyT' < T ande small enough, (e -1 (V — v),e~ 2 (p e — 1)) is bounded 
in L oo ([0 1 T r \,H a x F _1 (T d )). Finally when T = +oo, T e goes to infinity. 

Remark. The models that we consider are valid in domains without boundaries, and although 
stated in the space periodic case, we believe that our results hold true, with some technical adaptation, 
in the case of the whole space. 

2.2 Proof of Theorem I2TT1 
2.2.1 Heuristics. 

Let us introduce (v,p) the solution of the Euler incompressible system (|SJ) and corresponding pressure. 
Note that by taking the divergence of (JEJ) the pressure is given by the following: 

d 

Ap = divWjV 1 . 
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We will all along the paper use the following notation: for two vector fields u, v, 

d 

If v is solution to (EP e ), we introduce also 

v — v + evi, 
p = 1 + e 2 pi- 

We suppose for now that d = 2, and we take the curl and the divergence of the momentum equation, 
this yields 

d t cm\vi + v ■ Vcurlt>! = R 1 , 

d t divvi + v ■ Vdivt>i = — + R 2 , 

e 

„ _ divt>x 

d t pi + v ■ Vpi = h R 3 . 

e 

If we assume for now that v±, p\ and their spatial derivatives remain bounded (we do not specify in 
what sense yet), R^i = 1,2,3 are bounded terms. Hence the vorticity curb; is not affected by the 
electric field, and the vector u = (curl^i, divi>i, p\) evolves through 

d t u + v ■ Vu = — x u + R. 

e 

with R bounded, and k = (1, 0, 0). Under this form, the system looks like a rapidly rotating fluid (up 
to the remainder R, and the fact that t)^u) and the singular term - x u induces time oscillations 
of frequency e _1 , but does not increase the energy of the perturbation, allowing energy estimates, as 
long as the remainder R is under control. 

2.2.2 Reformulation of the system with new unknowns 

For a vector u G M, d , we denote u 1 , ...,u d its components. We define the new unknowns cu 1 , pi as 
follows: 

V-u = e/3i, 

p=l + e 2 pi = 1 + eA0, 
curl v = uj = u + ecui, 

with (v,p) the solution of (E), and u = curl v. 
We will use the following observation: 

Lemma 2.2 Let (p,v) be a solution to (EP e ). Then, the total momentum f p(t,x)v(t,x) dx does 
not depend on time. 
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Proof. The momentum equation can be rewritten 

dt(pv) + V ■ (pv <g> v) = pV0. 
We just have to show that the integral of the right hand side vanishes. For this, we use the identity 



(1 + e 2 A0)V0 = e 2 

that yields when integrated over T d . 
We note also 



V- (V0®V0)- ^V|V0| 2 



□ 



v — v + ev\, 

but v\ is not really an unknown since it can be obtained from the knowledge of o>i, fti, and e f pv 1 
J PqVq — J pv: Indeed when d = 2 we have 



difo + d 2 uji = Avl, 
d 2 /3i -diUi — Av \. 



(9) 
(10) 

In the 3 dimensional case we have equations (J9j fT0|) replaced by 

(V x V x v) + V(V ■ v) = Av, 

and thus 

(11) (V x uiY + di(3i = Av{. 

Note that when d = 2 the vorticity is scalar and when d ^ 2 it is a vector field of T d . We show now 
that v\ can be estimated from ui, and J pv. 

Lemma 2.3 Let (p, v) be the solution at time t of {EP € ) with initial datum (po, vo), let v be a solution 
at time t of (E) with initial datum vq. Then we have, for s G R ; 



v — v 



v — v\ 



< C \\V-v\ 



H"(T d ) 



+ II V X (v - v) 



\H s {T d ) j, 



v-v] 



< 



- illi^^dlf ||i2(Td) + ||V ■ v\\ H -i {Jd) + ||V x (v - v)||h-i(t*)) + J [PoVo-v ] j 

Proof. Let Ui be the unique vector field with zero average such that 

V x m = V x [v — v), 

V • m = V • v. 
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We have directly from identities (011011111) that ||mi||^+i < C^|| V • v ||#« + || V x (v — v)\\h*j , which 
is the first point of the lemma. The difference v — v — U\ = W\ is a constant vector field, and we have 



pv = J pv + J p(v — v) 

V+ (jj-l)v+ I (p - + Wi, 



which yields 



w\ = / [po^o -Vq]- / (p - l)v - Up - l)u 



and the result follows. 

□ 

We immediately deduce the following corollary: 
Corollary 2.4 Let (p,v) be a solution to (EP e ), let v be a solution to (E). Then, for any s > 0, 



\v-v\\h° < 



[p v - v ) 



+ c(i + i||h.)(||v ■ v\\ H s-i + ||v x (v-v)\\ H .-: 

+ C\\p - l\\ H s\\v\\ H s. 

Taking the curl of equation (jSJ) we recall the following identities: 

(12) d t u + (v ■ V)u + (V • v)u = when d = 2, 

(13) d t u) + (v ■ V)w + (V ■ v)u + (u ■ V)v = when d = 3. 

When d = 2 the (EP e ) system then becomes: 

(14) d t {u + ewi) + u • V(cj + ewi) = -(w + ewi)e/3, 

(15) ^e/3i + v • Vepi + 2ed i v j d j v\ + e 2 d i v{d j v\ = — - d^djv*, 

(16) d t e 2 Pl + i; ■ Ve 2 pi = -(1 + e 2 pi)eft, 



whereas when d — 3 one would have to replace equation (jT4j) by 

(17) 9f(a) + euj\) + t; • V(a> + eu;i) = — (u> + euj\)e(i — uj ■ Vv. 
Noticing that A<ft = epi, if we set 

(18) Pi = Pi - Ap, 
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we get the following system for d = 2: 

d t uji + v ■ Vo>i = — u/3i — eioifii — v\ ■ Va>, 



(19) 



Pi 



d t pi + v ■ Vpi = -— - e(pi + Ap)/3i - (d t Ap + v ■ VAp). 



2Vv x : - eVi;i : Vv ± , 



When d = 3, the first equation should be replaced by 

dfLUi + v ■ Vco>i = 

— vi ■ Vlj — uj ■ V^i — u; ■ Vv — eu)\j3\ — oo\ ■ Vt>i. 

2.2.3 Energy estimates 

We handle the energy estimates when d = 2 but the same result would hold when d = 3, just with 
more terms. The system can be written in the following way: 



(20) 
(21) 



dtu' + ^y<9;u e + R e u € = S e (u' 

i 

u e (0) = ug f 



where v is still the velocity, and where 






































1 

c 





The 'source' term S e is given by 

S e = 



—ujfli — euiPi — v\ ■ Vo; 
-2Vui : - eV^i : Vui 
-e(pi + Ap)/3i - (<9 t Ap + v ■ VAp) 



We apply <9 7 to equation (J27?j) . where 7 = (7 1 , 7 d ), and <9 7 stands for - f 7 q ^ d , with I7I = ^=1 7*. 
We get 



dxj ■■■9x 1 d 



where 



i=l |/i|>l,7>A« 



Then we have 
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Lemma 2.5 // |7| > d/2, then for S e and S e defined as above we have: 

ll^ e (^, -)IU2 (Td) < C(l + -)ll^i^i(T rf ) + e||u e (t, -)llLi.i ( Tr d )) 

and 

\\d-yS% -)|| L2(Td) < C(l + OHnhip*) + •)!^ W(Td) ). 

Proof. 

- Point 1: the proof is a straightforward adaptation of the proof of pQ p. 151. It is based on the 
following estimate: [Ij Proposition 2.1.2 p. 100: 

Proposition 2.6 Ifu,v E L°° D H s s E N, then for any 5,r), \5\ + \r]\ = s, one has 

WftuX&v)^ < c(||«MMI^ + NMMM- 

Applying this result to d ,x v l d n/ ~ IJ 'diU e , > 1, we obtain 

||<9 / V<9 7 ~ M <9 i u e || L 2 (Td) < C (||Vu t ||£»(Ti>i)||u e || £r M (T d) + llV^llHiTifrrfjIl^lliooppi)) . 
We know that || • IIl^cf*) < C\\ ■ \\g\f\nd) if \j\ > d/2. Thanks to IT0|) we have for any s 
(22) || Vvi\\ H s(jd^ <c(IM| H s (T d ) + ||/3l||^(Td)), 

thus || Vfi||^| 7 |( X d) < C'||u € || i3 -| 7 |( T d). Then using that v = v + ev 1 we conclude. 

- Point 2: We also know thanks to Proposition 12.61 that for s > d/2 

\\Vvt : VuiUff. < C||Vui|U«»||Vi;i||h'- 
It follows that for s > d/2, we have 

\\S e \\ H s <C(l+||u e ||^ + e||u e ||^) 
where C depends on the smoothness of the solution of ( J5j) . 

□ 

Having applied <9 7 to (J2*U|) . multiplied by <9 7 u e , and noticed that for any w, (w, R e w) = 0, we obtain 

^|d 7 uf + ^(^|d 7 uf) 
1=1 

= V ■ v ^|<9 7 uf + {d 7 S £ - £ e )<9 7 u e . 

Since ||V • f ||l°° = ||z,°° < Ce|| u<E ||/n if 7 > d/2, using Lemma 12*31 and integrating over T d we 
have, for any |7| < s, s > d/2: 

j t \\d 7 u% -)||i 2(Td) <c(l+ ||u e (t, Oll^j + e||u e (t, Oliver*)) ■ 
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By summing this over all multi-indexes |7| < s, we can conclude using a standard Gronwall's lemma, 
that if the solution v, p of Euler is smooth (see below for the smoothness required) on the time interval 
[0, T], for any X" < T there exists eo such that the sequence (u e ) e<eo is bounded in L°°([0, T'}, H S (R 2 )). 
Then we have 

u e = (curl^i, divf 1, pi — Ap), 
v = v + ev\ , 
p=l + e 2 Pl . 

We use Lemma 12.31 to get v from u, v, and from Corollary 12.41 the bound obtained on u implies 
a bound on e _1 (f — v),e~ 2 (p — 1). (Note that from the assumption on the initial data, we have 
\Jff Q v e -v \<Ce.) 

This proves Theorem 12.11 

□ 

Minimal regularity for the limiting field In order to perform our computations, we need to have 
at least f3\,p\,ui\ in L°°([0, T] xT d ) and thus to have an estimate on their norms in L°°([0, T], H s (T d )), 
with s > d/2. Therefore we need to control Then we need to apply <9 7 tofl9| with 

I7I = s > d/2, and we need to control \\d t Ap + v ■ VAp||# s (this is the 'worst' term). This implies to 
control ||u||fl-«+a. which requires v to be bounded in L°°([0, T], H s+2 (T d )) with s > d/2. If we take s 
integer, we ask s > [d/2] + 1, and v must be bounded in L°°([0, T], H^ d ^ +3 (T d )). 

Remark. Usually, modulated energy techniques only require a bound on || V^Hl^. Here we need 
one more derivative, since the 'div-curl' formulation of the system (performed in order to obtain 
energy estimates) is obtained by differentiation. 

2.3 Non-prepared initial data 

Here we obtain energy estimates in the general case of non-prepared initial data. What we mean by 
'general case' is the case of a generic smooth initial velocity, and smooth initial density, with finite 
energy, hence p — 1 will be of order e. 

We will see that the energy estimates are the same as in the case of prepared initial data, the 
asymptotic e — ► is then handled similarly to ^H], although the algebra in our case is quite simple. 
The solution will exhibit a good space regularity and a strongly oscillating behavior with respect to 
time. As explained in jTH] , the motion can be decomposed along a slow and a fast manifold; the slow 
manifold consists of divergence-free velocities with uniform density, and the fast manifold consists of 
potential velocities. Due to the rapid oscillations, the potential part of v will converge weakly to 0, 
and the divergence-free part will converge strongly to a smooth solution of the incompressible Euler 
equation. 

We still consider (p e , v e ) solution to (EP e ). For any vector field v, we introduce its solenoi'dal part 
and potential part, which is the pair (Uv, Vq) such that v = Uv + Vq with q periodic and V • Uv = 0. 
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A priori estimates We first rescale the density fluctuation 

P e = l + ep\. 

We use the unknown u e , given by 

(L) e = curlf e 
P\ 

(Note that we do not subtract Ap to p\ in this case.) 

We now restrict to the case d = 2, but one can check easily that the same results hold when 
d = 3. 

We write the equation followed by u e : 

(23) d t u e + rfdirf + R€ue = S e (u e ), 

i 

(24) u £ (0) = u$, 

where R e is a before, and the source term S e is now given by 

$■ = I -W : Vv e 

-PIP 6 



When s > d/2, proceeding as in Lemma f2. 51 we have the estimate 

II^IIh 8 < Cs|l u<: ||fp- 



Arguing as in the previous case, we conclude that, given s > [d/2] +2 and a sequence of initial 
data (Pq,Vq), with J p^ = 1, and such that 

e- x \\pl-l\\ Hs -i<C, 

there exists a sequence (p e ,v e ) of solutions to (EP e ) with initial data (pq,Vq), and (e _1 (p e — l),v e ) 
remains bounded in L°°([0, T], if s_1 x H s (T d )) for some T > 0, independent of e. 

Convergence We use the change of variable used for rotating fluids (see [21]), that removes the 
time oscillations: Considering the pair u e = (P e ,p\) such that 



G. Loeper, Euler-Poisson and Euler-Monge- Ampere 



13 



we have 

d t u £ + v e ■ Vu e = T e , 

with 

||^"Hl° o ([o,t],.h" 3 - 1 ) < C, 
|| l,e |U oo ([o,T],H s ) < C, 

||u e ||Loo([ T ] < C 

from the a priori bounds ( we still have s > [d/2] + 2). Hence, dtu e is uniformly bounded at least in 
L 2 , and we deduce classically that (f3 e ,p\) converges weakly to in [0,T] x T d . 

We assume now that TIvq, the soleno'idal part of Vq, converges weakly in L 2 to some limit v , hence 
it converges strongly in H s ' for s' < s. We check that d t uj £ is bounded uniformly on L 2 under our 
assumptions. Hence from the a priori bound, u e converges (if necessary passing to a subsequence, 
but see the remark below) in C([0, T],H S _1 (T d )) for all s' < s. 

Let v e be the unique vector field with zero average such that V x v e = u e . Then, Uv e , the 
soleno'idal part of v e is equal to v e + c e with c e a constant vector field. Since f p e v e = f PqVq, and 
since e~ 1 (p e — 1) is bounded in H s , we have limc e = lim J Vq — lim JU.Vq. Hence Tlv e = v e + c e 
converges in C([0, T],H s '(T d )) for all s' < s. 

Decomposing v e as v e = Vq e + Hv e , we have q e = /S.~ 1 f3 e that converges weakly to in L 2 ([0, T] x 
T d ). Hence, using the a priori bounds, in the vorticity equation 

dtuj 6 + v e ■ Vuj e = -f3 e u e , 

we can pass to the weak limit in L 2 ([0, T] x T d ) and state that u = lima; 6 satisfies in T>'([0, T] x T d ), 

d t uo + Uv ■ Vw = 0, 
u(Q) = uj , 

where is the limit of the soleno'idal part of v e . Moreover we have u = V x n^, hence Uv is a 
solution to the incompressible Euler equation, with initial data v the limit of the soleno'idal part of 
eg. 

Remark 1. From the regularity of vq, the solution (in the distribution's sense) to (E) with initial 
data nt> is unique. Therefore the whole sequence Uv e is converging. 

Remark 2. Here we did not introduce the solution of the limit equation (incompressible Euler in 
this case), and chose to argue by compactness. This method looks simpler, however, we obtain less 
informations concerning the 'rate of convergence' of the sequence (p e ,v e ) with respect to e. 

We gather those results in the following theorem: 

Theorem 2.7 Let s > [d/2] + 2, let (pq,Vq) be a sequence of initial data, such that J p € = 1. Assume 
that is bounded in H s (T d ) and that e _1 (Po — 1) is bounded in H s ~ 1 (T d ). Then 

1. There exists T > such that a sequence of solutions to {EP e ) with initial data {Pq,Vq) exists 
on [0,T] ; and the sequence (e _1 (p e - l),v e ) is bounded m L°°([0, T], H s ~ l x H s (Y d )). 
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2. The potential part ofv e , defined above, converges weakly to in L 2 ([0,T] x T d ). 

3. IfUvQ, the solenoi'dal part of Vq converges to some vq weakly in L 2 , then, for all s' < s, Hv e , 
the solenoi'dal part of v e , converges in C([0,T},H s '(T d )) to v the solution of (E) with initial 
data Vq. 



3 The Euler-Monge- Ampere system 

We consider here the following Euler-Monge- Ampere system denoted by (EMA e ): 

(25) d tP + V-(pv) = 0, 

(26) dtV + v-Vv= V ^~ X , 

(27) detD 2 tp = p. 

The last equation must be understood in the following weak sense: ip is the only (up to a constant) 
convex function with ip — \x\ 2 /2 being Z d periodic such that 

(28) VfeC°(T d ), [ f{V^)dp= I f(x)dx. 
This definition will be made precise in Theorem 13. II and Definition 13.21 



3.1 Geometric derivation of the Euler-Monge- Ampere system 

This derivation has been introduced in [S], [7j, [TH]. We reproduce it for sake of completeness, but 
the reader only interested in the proof of convergence may skip this section. 



3.1.1 The Euler equations of an incompressible perfect fluid 

The motion of an inviscid and incompressible fluid in a domain Q C M d is described by the Euler 
incompressible equation (E) that we recall here: 

d t v + (v ■ V)v = Vp, 
V ■ v = 0. 



Following Arnold (see |2]), we have a formal interpretation of the Euler incompressible equations: 
introducing G(Q) the group of all volume preserving diffeomorphisms of Q with jacobian determinant 
equal to 1, the Euler equations describe the geodesies of G(fl) with length measured in the L 2 sense. 
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3.1.2 Approximate geodesies 

A general strategy to define approximate geodesies along a manifold M (in our case M = G(Q)) 
embedded in a Hilbert space H (here H = L 2 (Q, M. d )) is to introduce a penalty parameter e > and 
the following unconstrained dynamical system in H 

(29) d tt X + ^Vx (d 2 (X, M))) = 0. 

In this equation, the unknown t — > X(t) is a curve in H, d(X, M) is the distance (in H) of X to the 
manifold M, i.e. in our case when M = G(Q), 

(30) d(X,G(Q))= inf \\X-g\\ H , 

g&G(n) 

finally, Vx denotes the gradient operator in H. This penalty approach has been introduced for the 
Euler equations by Brenier in It is similar-but not identical- to Ebin's slightly compressible flow 
theory ^3], and is a natural extension of the theory of constrained finite dimensional mechanical 
systems [22] • Actually if G(Q) were a smooth manifold, the result would be exactly the one of fT3~] . 
Theorem 2.7, but this is not the case, here because the L 2 metric is too weak. The penalized system 
is formally hamiltonian in variables (X, dtX) with hamiltonian (or energy) given by: 

e = l\\d t xf H + ±d 2 (X,G(n)). 

Multiplying equation (|2*9"j) by d t X, we get immediately that the energy is formally conserved. There- 
fore it is plausible that the map X(t) will remain close to G(Q) if it is close at t — 0. A formal 
computation shows that, given a point X for which there is a unique closest point tcx to X in the H 
closure of G(Q), we have: 

(31) V x (d(X, GO) = -^£){X - ir x ). 

Thus the equation (J29|) formally becomes: 

d a X + \(X-rr x ) = 0. 
e z 

To understand why solutions to such a system may approach geodesies along G(O) as e goes to 0, 
just recall that, in the simple framework of a surface S embedded in the 3 dimensional Euclidean 
space, a geodesic t — > s(t) along S is characterized by the fact that for every t, the plane defined by 
{s(i),s(£)} is orthogonal to S. In our case, d tt X(t) is orthogonal to G(Q) thanks to (}30|) and X(t) 
remains close to G(Q). 

Notation Since we intend to work on the flat torus T d we might consider Z d additive mappings, 
i.e. mappings such that 

\/peZ d , x(.+p) = x(.)+p, 



G. Loeper, Euler-Poisson and Euler-Monge- Ampere 



16 



as well as periodic mappings (i.e. mappings from T d into itself). 

Then given m an additive mapping,we denote by rh the naturally associated mapping on T d . The 
following polar factorization Theorem is a periodic version of jHj, it has been discovered independently 
by j201 and [32]. 

Theorem 3.1 Let X : ~R d — > M. d be Z d additive and assume that px = X#dx has a density in 
L 1 ([0, l] d ), then there exits an a.e. unique pair (V '(px^x) satisfying 

with (px & convex junction such that 4>x{ x ) — \ x \ 2 /2 is ^ d periodic, and ttx : M. d —>■ M. d additive such 
that Tlx is measure preserving in T d . Moreover we have 

\\X - 7Tx||L 2 ([o,i] d ) = 11^ ~~ ^xIIl 2 ^) = d(X, G(T 1 )) 
and if ipx is the Legendre transform of (f>x then 

IT x = VV>x ° X. 

Remark 1. The pair (4>x,i>x) is uniquely defined by the density px = X#dx. 

Remark 2. Important properties of the optimal potential: The periodicity of <px{x) — \ x \ 2 /2 
implies that V0x and Vipx are Z rf additive, and that ipx ~ \ x \ 2 /2 is a l so ^ d periodic. This allows 
the following definition: 

Definition 3.2 Let p be a probability measure on T d , then we denote <p[p] (resp. ip[p\) the unique 
up to a constant convex function such that 

'Pip] ~ I ' | 2 /2 is Z d periodic , 

VfeC°(T d ), [ f(y<j>\p](x))dx= f f(x)d P (x) 

Jl d JT d 

(resp. its Legendre transform). 

Remark. We recover thus that ?/>[p],0[p] will be generalized solutions of the following Monge- 
Ampere equations 

det D 2 ip = p, 
p(V0) det£> 2 = 1. 

3.1.3 Result 

We are now ready to state the main result of this section. 

Theorem 3.3 Let s e N with s > [d/2] + 2. Let v be a divergence- free vector field on T d , let 
(v,p) be a smooth solution of the Euler incompressible system on [0, T] x T d , with initial data 
vo and satisfying v G L°°([0,T], H s+1 (T d )). Let (v$, p e Q ) be a sequence of initial data, with f p e = 1, 
and such that (e -1 ^ — vo), e~ 2 (p € — 1)) is bounded in H s x H s ~ 1 (T d ). Then there exists a sequence 
(V,p e ) of solutions to (EMA e ) with initial data (v € , p e ) belonging to L°°([0, T e ], H s x H'' 1 ^)) with 
liminf e ^o T e > T. Moreover for any T' < T and e small enough, (e~ 1 (f <E — v), e~ 2 (p e — 1)) is bounded 
in L°°{[0 : T'],H s x H 3 ' 1 ^)). Finally when T = +oo, T € goes to infinity. 
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Before entering the proof of this result, we need some preliminary results concerning the lineariza- 
tion of the Monge- Ampere operator. 

3.2 Linearization of the Monge- Ampere operator in H s norm 

This section is devoted to the proof of the following Theorem: 
Theorem 3.4 Let p be a probability measure on T d , d < 3, let ip satisfy 

det D 2 ip = p 

in the sense of Definition ^. S\ Then, there exists eo > such that if \\p— ^H 2 {T d ) < eo, for any s£N 
with s > d/2 there exists C(s) that satisfies 

(32) \\D 2 ^ - I\\ H .(p) < C(s) ||p - 1 Hh'CT*)* 

(33) ||(A-0 -d)-(p- l)\\ Ham < C(s)\\(p - l)\\ 2 H s {Tdy 

3.2.1 Preliminary results 

We first state the following result obtained from jH] on the regularity of solutions to Monge-Ampere 
equation, adapted to the periodic case. This result will be the starting point of the proof of Theorem 

Theorem 3.5 Let p G C a (T d ) for some a > 0, with < m < p < M , be a probability measure on 
T d , let ip = ip[p] in the sense of Definition \3. Si Then ip is a classical solution of 

det D 2 ip = p 

and satisfies for any a' < a: 

(34) ||VV»(x) - x\\ L °o < C(d) = Vd/2, 

(35) \\D 2 i>\\ ca > < K{m, M, \\p\\ Ca , a, a'). 

Then we state a classical result of elliptic regularity that we will need during the course of the proof. 
It can be found in |14j . Theorem 9.11. 

Theorem 3.6 Let Vt be an open set in ~R d , u G Wf^iVt) D L P (Q), 1 < p < oo, be a strong solution of 
the equation 

d 

in Q where the coefficients satisfy 

a lj G C (O),f E L P (Q); 

m 2 < < Aier g R d , 

for i,j = l..d, with < A, A < oo. Then for any Q' CC Q, 

\\u\\w 2 'P(n') < C(||w|Up(n') + ||/IUp(n'))' 
where C depends on d,p, A, A, Q', Q and the moduli of continuity of the coefficients a*- 7 on Q' . 
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3.2.2 Proof of Theorem HOI 

Sketch of the proof We assume here d = 3. We recall that if) satisfies det D 2 ip = p. We first have 
to prove that p G H s implies D 2 ip G H s . We will proceed by induction. We recall first that for A, B 
two d x d matrices, we have the following expansion: 

det (A + tB) = det A + t trace (A^B) + o(t), 

where A com is the matrix whose elements are the minors of A, or co-matrix of A. Hence the elements 
of A com are polynomials of degree d — 1 in the elements of A. When A is invertible, we have 
A com = det A [A]- 1 . 

Differentiating s times the Monge-Ampere equation, and denoting M the co-matrix of D 2 ip, we 
will have 

trace (MD 2 d s if)) + T = d s p, 

where the first term contains the highest derivatives, and T will consist of products involving three 
derivatives of if). The order of each derivative will smaller or equal to s — 1, and the sum of the three 
orders will be equal to s. By a careful analysis, this product will be controlled in L 2 by ||Z} 2, ?/>||#s-i 
and HpII/p, using Sobolev injections. From Theorem 13.51 assuming a minimal regularity for p (i.e. 
the bound (|37j)). D 2 ip, and therefore M will be continuous elliptic matrices. Hence d s if> solves an 
elliptic problem, with continuous coefficients, and we will use the Theorem 13. 6l to obtain D 2 d s ip G L 2 . 

This intermediate step will be done in Lemma f3.7j using this a priori estimate and the continuity 
method, we will obtain the estimate (J32|) . 

Then, the expansion det(J + D 2 (p) = 1 + Atp + P(dij(p), where P is a polynomial in dij(p whose 
terms are of degree two or three (when d — 3), will yield (|3*3*jl . 

Rigorous proof We recall that if> satisfies 

(36) det D 2 if) = p, 

(37) ||p - 1|| H 2 < e , 

for some e to be chosen later. We suppose d = 3 and the proof can be reproduced in the case d — 2 
with minor modifications. The parameter eo is chosen such that (|3*7j) implies 

Ai < p < A 2 

for some Ai > 0, A 2 > 0. Note also that thanks to (|37jh p is in C a for a — |. Then from Theorem 
13.51 D 2 ip G C a ' with a' < a. Note also that since p G [Ai, A 2 ] and using equation (jSHJ), D 2 ip G C a ' 
implies that [D 2 ^] -1 G C a , and thus M lJ the co-matrix of D 2 ip is uniformly elliptic and in C a . 

We first prove by induction that if 7 G N d then p G i/' 7 ' implies D 2 if> G i?' 7 '. It can be checked 
during the proof that this bound will be uniform under the condition (|3"?j) for eo small enough. 

Lemma 3.7 Under assumption \3T\) , for any 7 eN d , pE implies that d^D 2 ip G L 2 . If moreover 
p G W^ e then d~ l D 2 if) G L 6 . 
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Proof. This lemma will be proved by induction. We first deal with the cases |7| = 0, 1, 2. 

The case 7 = is a consequence of Theorem 13.51 

For I7I = 1 we differentiate (jHHJ) with respect to x„, to obtain 

(38) M^d lJ (d^)=d u p, 

with the co-matrix of dijip. Then if d u p £ L 2 , by Theorem 13 .6\ d v ip £ W 2 ' 2 . If d v p £ L 6 we also 
get that d u ip £ W 2fi . 

For I7I = 2 differentiating once more with respect to x@ we obtain 

(39) M^d i5 {d v ^) + (dpW^dijid^), = d vP p 

still with M ij the co-matrix of <%0- Suppose that p £ H 2 , then W 2 ' 2 C W 1 '^ = W lfi if d = 3, 
and d v ip £ W 2,G . The term dpM l i is a sum of terms of the form dij{dpip)dkiip and the second term 
of the left hand side of (|39j) is thus bounded in L 2 . Then once again by Theorem 13.61 one gets that 
D 2 d uf3 ^ £ L 2 if d vP p £ L 2 . 

Moreover if D 2 p £ L 6 then d v p £ C a for some a > 0. Using (j38j) and Schauder interior estimates 
(see |Hj, Theorem 6.2.), we obtain D 2 d u ip £ C a ' . Thus (d /3 M ij )d i j(d u if;) £ From §^ and 

Theorem 13.61 we obtain d v pD 2 ip £ L 6 . 

As we just saw, Lemma T3. 71 is true for I7I = 0, 1, 2. We assume that it holds for all 7 with I7I < n 
for some n > 2. Take now |7| = n + 1 > 3, p £ il' 7 ', and apply <9 7 to (j36j) : 

(40) M ij dijd^ + * dijd^iljd kl d^tljd mp d r ^ = «9 7 p 

71+72+73=7 

I -y I — 1 > I Ti I > I T2 1 > I T3 1 

with * some constant coefficients. We call T the second term of the left hand side of (}4*0|) . Since 
peHW,pe W^' 1 ' 6 , and we have d a D 2 ip £ L 6 (T d ) for any \oc\ < n using the induction hypothesis. 
Therefore T £ L 2 and since d 1 p £ L 2 we obtain M^dijd^ip £ L 2 . Using Theorem 13.61 it follows that 
d^'D 2 ^ £ L 2 . 

Remember that I7I > 3 thus |7 3 | < ||7| < 7 — 2, and |7 2 | < \\^\ < 7 — 2. Since d = 3, we 
have i7 2 C C a for some a > and thus dkid^ip, d m pd 13 ip are in C a , moreover H l C L 6 and since 
|7i| < I7I — 1, dijd llr ip is in L 6 . Therefore T is in L 6 . If <9 7 p £ L 6 we have d^D 2 ^ in L 6 . Hence the 
lemma holds for all |7| < n + 1. This achieves the proof of Lemma f3. 71 

□ 

Now by induction on I7I we prove (j32j) and (jHHJ). From (j37J) we have ||p — 1||^2 < eo small. Take 
ip = \x\ 2 /2 + (p solution of det D 2 ip = p with <p periodic and J Td <p = 0. We begin to show that 
||(p||c-2, a for some a > 0, is controlled by ||p — 1||h 2 - Indeed, the periodic solution of 



det(/ + L>V) = p 
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can be built by the continuity method (see [Ej). Starting from p — l,(p — 0, we use the implicit 
function Theorem to obtain the solution cp t of 

(41) det(I + D 2 tp t ) = tp+ (1 -t). 
For this we differentiate (}4~Tj) with respect to t, to obtain 

(42) Mi j d ij a t tp t = p-1, 

for t G [0, 1], where M t is the co-matrix of I + D 2 ip t . We know, from the a priori estimate of Theorem 
I3.5[ that for p G C 1 / 2 , D 2 cp t and therefore M l t 3 are C a elliptic matrices for all a < 1/2. 

To see why f|42|) indeed admits a unique (up to a constant) periodic solution, we recall that M is 
the comatrix of a Hessian matrix, therefore it is 'divergence- free': 

V* G [l..d], ^M ij = 0. 

j 

Hence equation (|4*2*j) can be rewritten in divergence form 

^(M^^) = p-1, 

and the operator C = M\ 3 diy is a self adjoint operator on if 1 (T d ) and induces a bounded coercive 
bilinear form on H\^f. d \ where the subscript means that we impose the mean value to be 0. Then 
the existence/uniqueness of a solution to f|42H in follows by Lax-Milgram Theorem. 

Hence, d t (pt is the unique (up to a constant) periodic solution of the above elliptic problem, and 
from Schauder interior estimates we obtain 

||<9^I|C 2 .« < C\\p— l|| C l/2 

< Cllp-lHaa, 

uniformly in t G [0,1], and finally for a < 1/2, 

ll^llc 2 ' — c II p _ i||h 2 - 

Then we have 

det(J + D 2 tp) = 1 + Aip + Rijdijip 

where R is a symmetric matrix whose coefficients are polynomials in dijip of degree larger or equal 
to 1. The norms ||i2jj||ca are controlled by ||p — l\\c<* < to, hence, for eo small enough, the matrix 
6ij + Rij is uniformly bounded, elliptic, and C a continuous. Since satisfies 



(43) 



(Sij + R^dijip = p - 1 
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it follows from Theorem 13 . 61 that Hc^^H^ < C\\p — 1||^2 and this proves (}3*2*|) for 7 = 0. 

If I7I = 1, we have 

(M ij )d i:j d u <p = d uP 

with M uniformly bounded, elliptic and C a continuous. For the same reasons we have ll^c^ll^ < 
C\\d v p\\ L 2. 

If I7I = 2, we do as in the Proof of Lemma l3~Yl using (J38| l39|) and keeping track of the bounds, 
we get 

\\dijd v pip\\ L i < C||9^p|| L 2 + C\\p- 1\\ 2 H 2 

< C\\ P -1\\ H 2, 

with C uniform under the assumption of Lemma 13.71 

If |t| > 3, we go back to equation (jHH): T is a sum of terms which contain all a product of at 
least two derivatives of ip of degree higher than 3. Since D 3 ip = D 3 (p, we have 

dijd^ijjdkid^ipdmpd^ip = dijd^ipdkid^ipdmpd 13 ^. 

We assume by induction that II^-^H^m-i < C\\p — l||#| 7 |-i. Since p — 1 G if' 7 ' we also have that 
D 2 ip G if' 7 ', with a uniform bound thanks to Lemma f3. 71 We remember that when I7I > 3, we have 
I72I < |t| ~ 2, |7 3 | < I7I — 2, therefore, using the injection of H 2 in C 1//2 when d = 3, dkid 12 Lp and 
d mp d ri ijj are uniformly bounded in L°°. We obtain that 

Wd^^dMd^d^d^W^ < C\\p- l\\ HH -i. 

Then <9 7 y? satisfies 

(M ij )d i3 d^ = d^p-T 

thus, using Theorem 13. 6[ 

H^DVIb < CdlP-lllffWM + ll^plUO 

< c(llp-lll H l T l), 

and we conclude that 

\\I?iI>-I\\h'<C{s)\\p-1\\h. 
for s G N, s > 2 and under condition (|3*7jl: thus (|3~2|l is obtained. 

Using Proposition 12.61 and the fact that dmpd 13 ^ is uniformly bounded in L°° thanks to Lemma 
13.71 we can also obtain that 

< C\\p-l\\ 2 Hhl 
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for I7I > 3. (When I7I = 2, the estimate holds also, but not using Proposition 12.61 ) Therefore, for 
all 7, we have ||T||£2 < C\\p — 1||^| 7 |. 

To conclude (}3*3*j) . we now write det(J + D 2 ip) = p under the form 

A<p = p-1 + P(d ij (p), 

with P consisting of products of two or three second derivatives of (p. Hence, under assumption ([37)1 . 
using Proposition I2.6| we have for s > d/2, \\P{dijip)\\H>> < C||Z? 2 <^||| fs . Using the bound (J32|) . we 
conclude that for s G N, s > d/2, 

\\A<p-(p-l)\\ Hs < \\p-l\\ 2 Hs . 

and Theorem 13.41 is proved. 

□ 



3.3 Energy estimates and proof of the convergence 

The proof of the energy estimates for Euler-Monge- Ampere is much inspired from the proof of The- 
orem 123 for the following reason: by taking the divergence of equation (}26|) one gets: 

dt(y -v) + v v(v • v) + d^djv* = ^ ~ d . 

Suppose that p is close to 1 at an order e 2 as is the case for Euler-Poisson, we guess (from Theorem 
13. 4j) that we have the following: 

ij> = \x\ 2 /2 + e 2 cp, 

p = det D 2 ip = 1 + e 2 A<p + 0(e 4 ), 

and thus = d + e 2 Aip = d + p — 1 + 0(e 4 ). Therefore we expect that 

d t (V ■ v) + v ■ V(V • v) + dyay = + 0(e 2 ), 



and that the technique of Theorem 12.11 will apply. 

Before performing the div-curl decomposition for the energy estimates, we need to establish the 
analog of Lemma \2. 21 in the present case, so that Lemma \2. 31 and its Corollary 12.41 hold. 



Lemma 3.8 Let (p,v) be a solution to (EMA e ). Then, the total momentum J p(t,x)v(t,x) dx does 
not depend on time. 



Proof. We proceed similarly as in the proof of Lemma 12.21 We need to establish that 

pit, x)(Vif)(t, x) — x)dx = 0. 



id 
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For this we use Definition 13.21 For f(x) = Vip(x) — x, we have 

/ Pf = [ /(V0) 

JT d JT d 

= [ V^(V0)-V0 
x — V0 

= o, 

where we have used at the third line that, for if), <f> Legendre transform of each other, Vip(V(f)) = id, 
and at the last line that — |x| 2 /2 is periodic. 

□ 

Hence we have shown that one can retrieve v from the initial value of J pv, and V • v , cur If . 




General framework We perform the same div-curl decomposition as in the Euler-Poisson case. 
We then express the difference between the solution of (EMA e ) and the limiting solution : either the 
solution of (E) or the solution of (EMA e ). After having applied a proper scaling to this difference, 
our solution is now described by a vector u whose first component (that can be a vector if d = 3) is 
the rescaled vorticity, and whose last two components are a rescaled divergence and rescaled density 
fluctuation. For this perturbation we will obtain 

d t u e + *Mu e + R e u £ = Q £ (u e ), 

i 

where we still use 

/ \ 
W = -± . 

V o \ o e ) 

For the source term Q e , we have < C(l + ||u e ||# s + <5 e ||u e ||#s) where 5 e goes to when e goes 
to 0, and the constant C depends on the regularity of the limiting field. 

Regularity of the limiting field The form of the source term will vary under the circumstances, but 
the general idea is that in order to bound Q e in H s , we will need the limiting velocity to be bounded 
in H s+2 and the limiting density to be bounded in H s+l . Remember that the H s norm of u controls 
the norm of (p, v) in H s x H s+1 , thus the limiting field must have one more derivative bounded than 
the order of the energy estimate. 

A Gronwall's lemma then yields a control on the perturbation that holds on a range of time [0, T e ], 
where T e — > T, T being the time of existence of a smooth solution for the limiting equation. 



G. Loeper, Euler-Poisson and Euler-Monge- Ampere 



24 



Convergence to Euler, two dimensional case Doing the same change of variables as in the 
proof of Theorem 12.11 



V - v = e/3i, 

p = 1 + e 2 pi, 

curlt> = a; = to + ec^i, 



we obtain: 

(44) d t (u + eui) + v ■ V(uj + eu x ) = -(u> + eu^efi, 



(45) dtefa. + v ■ Ve(3 x + 2ed i v j d j v\ + eS^'S,-*;} = — - - o^'fi^ 



(46) d t e 2 Pl + u • Ve 2 pi = -(1 + e 2 Pl )e^. 

Now we define E by 

A4>-d = e 2 pi + e 4 2, 

and from Theorem 13.41 inequality (J3*3*j) . we have, if s > 2, ||H||h« < The system can here 

be written in the following way: 

d t u e + J^y<9;u e + iTu £ = S e (u e ) + 
u £ (0) = u^, 

still with 



u 







r 







Pi 
















V° 


1 

c 





with the same S e as in the Euler-Poisson case, pi — pi — Ap, and with 



r = I eS 



We have ||V e ||H« < Ce(l + Hu^H^), for s large enough. Then the energy estimates are the same as 
in the first proof, the solution u e satisfying a control of the form: 

|K||^<c(i + KH^ + 6|His r .) 

and the same conclusion holds true. Then from Corollary 12.41 v — v, p — 1 can be retrieved from u e , 
and we obtain the expected conclusion. 

□ 
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Convergence to Euler, three dimensional case In the 3-d case, equation should be re- 
placed by 

d t (u) + eui) + v ■ V(a> + eui) = — (u) + ea>i)e/3 — uj ■ Vt>. 

Note that the vorticity equation is the same as in the Euler-Poisson case. This change would not 
affect the energy estimates. 

Higher order approximation Here we prove that the the Euler-Poisson system and the Euler- 
Monge-Ampere system are closer as e goes to than Euler-Poisson and Euler. We fix s > so = 
[d/2] + 2. For vq a H s+2 smooth divergence-free vector field, we consider (v,p) a solution of the Euler 
incompressible system (JBJ) such that v G L°°([0, T], H s+2 {T d )) for some T > 0. We consider also a 
sequence {v e ep1 p ep ) of solutions of the (EP e ) system with initial data (v e Q) p e ep ^) (with f p$ = 1) such 
that e-\v% p - v),e~ 2 (pl p - 1) is bounded in L°°([0, T'], H s+1 x H s (T d )), for any < V < T, if e 
is small enough. Thanks to Theorem 12 . 1| and from the regularity assumption made on v, such a 
sequence exists for any sequence of well prepared initial data. 

Theorem 3.9 Let s£N with s > [d/2] + 2. Let v,vl p ,p e ep be as above. Let (vq,pQ be a sequence of 
initial data such that (e _2 (t>o — vl p ), e~ 3 (p e — p £ ePiQ )) is bounded in H s x H s ~ 1 (T d ). Then there exists 
a sequence (v e ,p € ) of solutions to (EMA e ) with initial data (vQ,p e Q ), belonging to L°°([0, T e ], H s x 
H s ~ 1 (T d )) , with lim inf T e > T. Moreover for any T' < T and e small enough, the sequence 
(e- 2 (v e -v e ep ),e~ 3 (p e -p e ep )) %s bounded m L°°([0,T% H s x H*- 1 ^)). 

Remark. We see here that the difference between solutions of {EP e ) and (EMA e ) is of order e 3 
for the density and of order e 2 for the velocity whereas the difference between solutions of {EP e ) (or 
(EMA e )) and Euler was of order e 2 for the density and of order e for the velocity. 

Proof. We give the proof when d = 2, the proof would be the same when d = 3, just with more 
terms. We introduce (v ep , p ep = 1 + t 2 pi) solution to {EPA, and {(3 ep , u) ep ) = (V ■ v ep , V x v ep ). Then 
we set 

v = v ep + e 2 v 2 , 
V-v = (3 ep + e 2 (3 2 , 

P = Pep + e 3 p2, 

curk> = uj = oj ep + e 2 uJ2- 

The system {EMAA now reads: 

(47) d t {uj ep + e 2 u 2 ) + v ■ V{u ep + e 2 u 2 ) = -{u) ep + e 2 u 2 ){(3 ep + e 2 (3 2 ), 

(48) d t {(3 ep + e 2 (3 2 ) + v ■ V(/5 ep + e 2 (3 2 ) + V{v ep + t 2 v 2 ) : V{v ep + t 2 v 2 ) = 

(49) d t {p ep + e 3 p 2 ) + v ■ V{ Pep + e 3 p 2 ) = -{p ep + e 3 p 2 ){p ep + e 2 (3 2 ). 
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We still define 5 by 

- d = p - 1 + e 4 H, 

and from Theorem 13.41 we will have 

||"||^(T d ) < 4 \\p — l\\ 2 H s(jd) < C(l|Pl||^(T d ) + f2 II P2H 

(we use the notation p ep = 1 + e 2 pi). 
Setting 




u 



we obtain that 

d t u e + v Vu e + R e u e = T e , 
with R e as before and T t defined by 

-V 2 ■ Vw ep - (3 ep U 2 ~ (3 2 U ep ~ £ l U<l{h 

% v - 2V 



,~c 2 > 

T e =1 -v 2 - V/? ep - 2Vv ep : Vv 2 - e 2 Vv 2 :Vv 2 + E 



-ev 2 ■ Vpi - P ep p2 ~ efopi ~ £ 2 p2p2 
Using again Proposition 12.61 as in Lemma \2. 51 we obtain that, for s > d/2, 

\\Te\\H s (T d ) < Cs(l + || u<: ||_ff s (T d ) + e ll U<E ||/p(T d ))' 

where the constant C s is controlled by ||f e p||H s + 2 , (still with p\ = e~ 2 (p ep — 1)). From 

Theorem 12.11 these quantities are controlled for 0<t<T'<T,T being the time on which the 
solution of (E) is smooth. Hence we have by Gronwall's lemma a bound on ||u e ||£oo(r 0iT /] H sy 

Arguing as in the previous proofs, and using Corollary we obtain that (v 2 , p 2 ) remains bounded 
in L°°([0, T"\ , H s+1 x H s (T d )) for any T" < V and for e < e small enough. It follows that (e- 2 (v ep - 
v e ma),£~ 3 (Pep — Pema)) remains bounded in L°°([0, T"], H s x i7 s_1 (T d )). This achieves the proof of 
Theorem EUfl 

□ 



3.4 Non-prepared initial data 

In this case, we obtain exactly the same result as for Euler-Poisson, using the same techniques. We 
follow closely section 2.3, and we only have to estimate the additional source term that will appear 
in the equation followed by /3i, due to the Monge- Ampere coupling. We recall that det D 2 ip = p, and 
we will have to estimate the difference 

i [(Aip -d)-(p- 1)] 

in H s when we know that e _1 (p — 1) is bounded in H s . Thanks to Theorem 13.41 we conclude that for 
s > 2 this term is controlled by || e 1 (p — l)||# s = ||pi||h s - Hence the energy estimate can be handled 
similarly just with an additional term, and the conclusion remains true. 
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Theorem 3.10 The Theorem~2~7l holds also when replacing the (EP e ) system by the (EMA e ) system. 



□ 
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